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Abstract. Treating the motion of a dust particle suspended in a liquid as a random walk, Einstein in 1905 derived an equation 
describing the diffusion of the particle's probability distribution in configuration space. Fokker and Planck extended this work 
to describe the velocity distribution of the particles. Their equation and its solutions have been applied to many problems 
in nature starting with the motion of Brownian particles in a liquid. Chandrasekhar derived the Fokker-Planck equation for 
stars and showed that long-range gravitational encounters provide a drag force, dynamical friction, which is important in 
the evolution of star clusters and the formation of galaxies. In certain circumstances, Fokker-Planck evolution also describes 
the evolution of dark (invisible) matter in the universe. In the early universe, the thermal decoupling of weakly interacting 
massive particles from the plasma of relativistic leptons and photons is governed by Fokker-Planck evolution. The resulting 
dissipation imprints a minimum length scale for cosmic density fluctuations. Still later, these density fluctuations produce 
stochastic gravitational forces on the dark matter as it begins to cluster under gravity. The latter example provides an exact 
derivation of the Fokker-Planck equation without the usual assumption of Markovian dynamics. 



INTRODUCTION 

While working in 1905 at the Swiss Patent Office in Bern, Albert Einstein wrote two papers and his PhD thesis on 
the topic of Brownian motion id 0] • These papers and several later works are translated into English in Ref. [4]. 
Einstein's work on Brownian motion helped Jean Perrin confirm the existence of atoms and molecules [5] and laid the 
foundations for stochastic processes in statistical physics, a topic of great importance a century later. 

Brownian motion is the erratic movement of macroscopically small bodies suspended in a liquid. The phenomenon 
was described in 1827 by biologist Robert Brown |6|. Brown recognized that random motion occurred equally for 
living and nonliving particles. For six decades the phenomenon attracted little attention, although several authors 
suggested that the motion was due not to convective currents in the liquid but instead might be caused by molecular 
collisions. This hypothesis gained credence with the experiments of Gouy |7], which stimulated increased attention to 
the subject. Einstein — and independently Sutherland 1 8 ] — provided the first quantitative theory of Brownian motion, 
followed closely by Smoluchowski |9]. 

An excellent brief account of the research inspired by Brownian motion is given in Ref. lHoll . Pais's scientific 
biography of Einstein lUlll is the definitive source for the history and context of Einstein's work. The current article 
focuses on the role of Brownian motion in the development of kinetic theory and its application to weakly collisional 
systems in astrophysics. 

EINSTEIN'S ANALYSIS OF BROWNIAN MOTION 

Einstein recognized that a system of particles suspended in a liquid would have osmotic pressure and would undergo 
diffusion like a mixture of gases. As he did so often in his major works, he evaluated the dynamical equations at thermal 
equilibrium. Consider a suspension with uniform temperature having a number density gradient — for example, 
a teaspoon of sugar being dissolved in water. Treated as a fluid, the suspension obeys continuity and momentum 
equations. If the mean free time for collisions with the solvent atoms is sufficiently short, then the suspension will 
be in a quasi-equilibrium state even with nonzero density gradient and fluid velocity. One-dimensional equilibrium 
implies 



1 . Force balance: / — dp/dx = 

• Ideal gas law (van 't Hoff): p = nRT /Na (p = osmotic pressure, n = number density of suspended particles, 
R = gas constant, Na = Avogadro's number) 

• Stokes drag force: f/n = 6nriav (rj = viscosity, a = particle radius, v = velocity) 

2. Number flux balance: nv — D x dn /dx = (D x = diffusivity) 

The pressure force and diffusive flux are both proportional to dn / dx, which may be eliminated to give the Einstein- 
Sutherland relation, 

D x = -^—. (1) 
6nr\aNA 

If D x , rj, and a can be measured, this relation allows the determination of Avogadro's number. Modern physicists are 
so used to the ideal gas law p = nkT that we forget the atomic origin of Boltzmann's constant, k = R/Na- 
In the limit of negligible advective velocity v, the number density obeys the diffusion equation 

dn d 2 n 

whose Green's function solution in an unbounded domain is 

n(x.t) — expf— — — ] . (3) 

We see that diffusion leads to a steadily growing mean square displacement for each Brownian particle, 

(x 2 ) = 2D x t . (4) 

This is the same law as a random walk because diffusion is equivalent to an ensemble of random walks. A Java applet il- 
lustrating Brownian motion as a random walk is available online at |http : / / www . aip . or g/ hi s t ory / e ins te in /brownian . hi 

These results allow one to measure Avogadro's number Na by measuring the temperature and viscosity of the solute 
and the size of the Brownian particles, and then measuring the diffusivity from (x 2 )(t) for Brownian motion |12]. 
Perrin achieved the technical breakthrough of uniform-sized Brownian particles by repeated centrifugation of resin 
particles over a period of months. His work on sedimentation and the "discontinuous structure of matter" (i.e. the 
existence of atoms, in the stilted English of the Swedish Academy) led to the Nobel Prize in Physics in 1926. 

Einstein was unable to measure the size of Brownian particles this way. Instead, in his PhD thesis, Einstein cleverly 
deduced how the viscosity of a suspension varies with the volume fraction of the suspended particles, which is 
proportional to ct'Na- If the diffusivity is known, then equation Q leads to an estimate of Na- Einstein considered 
a sugar-water solution, boldly treating the sugar molecules themselves as Brownian particles, even though their 
individual motions were invisible. Measurements of the diffusivity and viscosity as functions of concentration led 
him to an estimate of the sucrose molecule size and Avogadro's number: a = 6.2 A, Na = 3.3 x 10 23 . The calculation 
contained an arithmetic error in the viscosity calculation, which was pointed out several years later by a student of 
Perrin's. When Einstein corrected the error |2J], he obtained a — 4.9 A, leading to the much more accurate result 
N A = 6.6 x 10 23 . 



KINETIC THEORY OF BROWNIAN MOTION 

Brownian motion can be described in two complementary ways: random walks and diffusion. The former approach 
follows individual particle trajectories while the second follows a distribution function. In this section we begin with 
the Einstein-Smoluchowski theory of random walks and conclude by exploring diffusion in velocity space with the 
Fokker-Planck equation. 

The simplest description of Brownian motion is a sequence of impulses (instantaneous velocity changes) separated 
by ballistic motion. After steps of duration At = t /N starting from position xo, 

Af 

X = X + V ! Af ' 



Suppose that the velocities have zero mean and are statistically independent (a Markov process), with 



(v,-)=0, (v,^} = o?Ify 

where I is the unit tensor. Then by the Central Limit Theorem, as N = t/(At) — ► °°, x — Xq becomes Gaussian with 
covariance 

((x-Xo)(x-Xo)) =N(a v Atfl = 2D x tI , D x = l -a 2 {At) , (5) 

recovering the result of equation ®. 

Random walks in position are unrealistic because the velocity cannot change instantaneously. Assuming a random 
walk in velocity with zero-mean, statistically independent accelerations a,-, repeating the above derivation gives 

of = 2D v t, D Y = l -a 2 a {At) . (6) 

In thermal equilibrium, a 2 = kT /M where M is the mass of the Brownian particle, yielding the unphysical result 
kT = 2MD v t — ► °° as t — ► ». It was invalid to assume zero mean acceleration. As a Brownian particle's velocity v 
increases, it sees a larger flux of background particles moving opposite its velocity, and collisions with them transfer 
a net momentum proportional to —v. (Stokes drag is the macroscopic equivalent.) Assuming a linear drag force and 
discrete Markovian dynamics, averaging over the acceleration for a given velocity gives 

2D 

(a i ) = -rv,- ) ((a,.-(a I -))(a ; -{a ; -))) = ^I%. (7) 

Now, v, + i = v, +a,(Af). Averaging over velocities, in equilibrium we must have (v,) = (v,+i) = and (v,v,) = 
(v!+iv,+i) = {kT /M)\. Working to first order in At, and imposing the Markovian condition ((a; — (a,))v,) = 0, we 
obtain the important result 

This formula is an example of the fluctuation-dissipation theorem: in equilibrium, diffusivity (describing fluctuations) 
is proportional to damping. At first glance, the Einstein-Sutherland relation Q seems to violate this because the 
diffusivity is inversely proportional to the viscosity. However, the Einstein-Sutherland derivation is valid only in the 
limit of overdamped motion, for which D v = J 2 D X with y = 6nrja/M. 

A continuous description of Brownian motion in velocity space was made possible by the introduction of a stochastic 
differential equation by Langevin in 1908 IU3I1 : 

d\ 

- = -yv+r(t), (9) 
where T is a stochastic force with mean and covariance 

(T(0)=o, (r(t)r(t-t')) = 2D v 8(t-t')i. (io) 

Langevin multiplied equation (|9} by x and averaged to derive the Einstein-Sutherland relation, which he showed is 
valid only for yt ^> 1 . However, the Langevin equation has a much wider applicability to stochastic phenomena in 
physics and other disciplines. 

The Langevin equation focuses attention on particle trajectories. It is often more convenient to describe a system 
statistically using distribution functions. A lucid translation between these two descriptions was provided by Klimon- 
tovich 1 14] and Dupree A system of particles has one-particle velocity distribution function 

/(v,O=E0(v-v,(O), <Hv) = S 3 (v). (11) 

i 

Here, i labels the particles and the Dirac delta function is a unit-normalized distribution. Evolving the system for a 
time At and Taylor expanding this distribution gives 



dt 



d<j) d\i At d\i d d\i 
d\ dt 2 dt d\d\ dt 



-0(At) 2 . (12) 



Now taking the ensemble average and using equations (|9jl and dlOt yields the Fokker-Planck equation fl^l : 

df d f df\ ykT 

-^- = —- yv/ + D— , H=- — I. (13) 
dt dv V dv J ' M 

This important equation describes diffusion (heating) and drag in velocity space. The fluctuation-dissipation theorem 
— the linear relation between diffusivity and drag — is necessary to ensure the correct equilibrium solution 

M \ 3/2 ( Mv 2 \ 



f W = n [2Mr) ^{-ikf)- (14) 

When the distribution function depends on position as well as velocity, the Fokker-Planck equation generalizes to 
the Kramers equation fl7ll (often called Fokker-Planck): 

°f- d f+\ d f +g d f - d (a/ yvf D d f) (15) 
dt dt dx d\ d\ \ d\ J 

Here g is the fluid acceleration, while the quantities A (drift), y (drag), and D (diffusivity) are called transport 
coefficients. Often, but not always, they are independent of v. (The distinction between g and A is then purely 
conventional.) Note that Df/dt represents the transport (advection) of particles in the (x,v) phase space along 
characteristics dx/dt = v, d\/dt = g. The quantity J = A/ — yvf — D • df/d\ is a flux density in velocity space. 
The total particle number f fd^xd^v is conserved. 

Equation dl5> and its relatives are called advection-diffusion equations. They have widespread applications in 
plasma physics, astrophysics, and other disciplines fill fTotl . Although our analysis began with Brownian motion, 
the advection-diffusion equation can describe many systems in which particle interactions play a dynamical role, 
including weakly collisional gases. The advection-diffusion equation and its transport coefficients must be derived for 
each application from more fundamental dynamics, or justified phenomenologically. 

The rest of this article will discuss applications of the advection-diffusion equation in astrophysics. 



BROWNIAN MOTION OF STARS 

Our galaxy contains more than 100 globular clusters, dense balls of 10 5 or more stars formed early in the galaxy's 
history. Some of these clusters contain many more X-ray sources than are found among a comparable number of 
stars in low-density environments elsewhere in the galaxy 1 20] . The X-rays arise from gas transferred from a close 
companion and accreting onto a compact object (neutron star, white dwarf, or black hole). The implication is that some 
globular clusters have many more close binaries (per unit mass of stars) than the rest of the galaxy. Why? 

A plausible answer to this question was suggested more than 30 years ago l2lll . Gravitational scattering between a 
binary and a third star can transfer energy to the third star, resulting in a more tightly bound (hence more compact) 
binary. If the internal velocity of the binary is greater than the typical speed of other stars (i.e., it is a "hard" binary), 
encounters preferentially remove energy from the binary (causing it to "harden"). Once a stellar binary is sufficiently 
hard, stellar evolution can form a compact object which accretes from its companion. 

When energy is removed from a binary the relative velocities of the two stars increases. This is true for any self- 
gravitating system close to dynamical equilibrium, as can be seen from the classical virial theorem. Let K and V be the 
total kinetic and potential energy per unit mass of the cluster, respectively. The virial theorem for an inverse square law 
force states (2K + V) ~ where the angle brackets denote a time average over a timescale longer than the characteristic 
time for stars to cross the cluster. Using the virial theorem we can evaluate the temperature-dependence of the specific 
heat. The total energy per unit mass is E = K + V and the "temperature" is proportional to the kinetic energy per unit 
mass. 1 The specific heat is then 

dE d , . dK 

_«_(*+y) = (16) 

where the virial theorem has been used. Self-gravitating systems have a negative specific heat and are therefore 
thermodynamically unstable [22]. This is true for any central force with two-body potential proportional to r" with 
—2 < n < 0. It is also true for black holes in general relativity (for which the temperature is the Hawking temperature.) 



The velocity distribution for a system of bodies with only gravitational forces is generally non-Maxwellian. 
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FIGURE 1. Spherically averaged density profile of a N = 1000 star cluster (the average of 5 clusters) evolved starting from a 
Plummer sphere with core radius r = 3tt/ 16. (The total mass and energy of the cluster are 1 and —1/4, respectively; the crossing 
time is f cl - = 2y2 and the relaxation time is 1 K \ = 0.1Nt CT /lnN = 41.) Two-body scattering transfers heat out from the core of the 
cluster causing it to shrink and become increasingly dense until hard binaries form. 



The thermodynamic instability operates on the timescale for particle interactions to thermalize the cluster. Suppose 
that the core of a cluster shrinks slightly. As a consequence of the virial theorem the core heats up. Scattering by 
gravitational interactions conducts heat outward, causing the core to lose energy and thereby contract further. The 
outer part of the cluster expands with the addition of energy. Heat conduction in a stellar dynamical system is a 
diffusive process governed approximately by the Fokker-Planck equation 112 311 . 

Qualitatively, the kinetic theory for a stellar dynamical system describes fluctuations caused by two-body scattering 
superposed on mean-field dynamics. The rate for two-body scattering can be estimated using the Rutherford (Coulomb) 
cross section for scattering of bodies with relative speed v: 



where we have assumed equal mass stars; the "Coulomb logarithm" In A accounts for long-range interactions, which 
dominate the total cross section for momentum transfer. The drag coefficient y °< v~ 3 for gravitational interactions; for 
a nearly Maxwellian distribution with velocity dispersion a v the velocity-space diffusivity D w ya; 2 . For the densest 
globular clusters, y l ~ 10 9 yr is much less than the age of the clusters, implying that two-body scattering has had 
ample time to drive the collapse of the core and, plausibly, the formation of compact binaries which then become 
X-ray sources. 

Modern computing power makes it almost possible (with approximate treatments of stellar evolution and tidal 
interactions) to study globular cluster dynamics directly using full numerical simulation 1 24] . A more modest example 
illustrating the gravothermal instability is shown in Figure^using exact pairwise forces and a fourth-order symplectic 
integrator 1 12511 . 



Most of the mass in the universe is in some form other than atomic ("baryonic") matter. This matter is invisible — it 
neither scatters nor absorbs detectable amounts of electromagnetic radiation — and has been dubbed "dark matter." 
Dark matter has been detected only through its gravitational effect on atomic matter and light. Its existence has been 
suspected for 70 years and known with confidence for about 25 years. However, the exact nature of dark matter is still 
unknown. 




(17) 



DARK MATTER AND BROWNIAN MOTION IN THE EARLY UNIVERSE 



The most natural hypothesis is that dark matter consists of particles or fields not yet detected in the laboratory. As- 
trophysics and cosmology provide strong constraints on the nature of this substance. The particles (or field excitations) 
must have sufficiently small thermal speed to provide the seeds for galaxy formation. This rules out standard model 
neutrinos. It cannot be atomic matter in any form without having overproduced helium during big bang nucleosyn- 
thesis. While proposals have been made that dark matter is an illusion arising from modified gravity 1 26], the only 
known theoretically and experimentally consistent implementation of modified gravity without dark matter l27ll is far 
less economical and less well tested than the dark matter hypothesis. 

Dark matter is natural in extensions of the standard model |28], with two leading candidates: axions and WIMPs 
(weakly interacting massive particles). Axions are hypothetical ultra-low energy excitations of the 0-vacuum of QCD 
that arise naturally as a solution of the strong CP problem (i.e., why the strong interactions conserve CP, or why the 
neutron has a tiny electric dipole moment). They are spin-0 particles of mass 10~ 6 to 10~ 2 eV and despite their light 
mass have negligible thermal speeds because they form a Bose-Einstein condensate. 

WIMPs are popular dark matter candidates because supersymmetry and other extensions of the standard model 
of particle physics call for particles with approximately the correct mass and cross section to produce the observed 
abundance of dark matter 12911 - The favored dark matter candidate is the lightest neutralino %°, the supersymmetric 
partner of a linear combination of the photon, Z°, and Higgs bosons. The neutralino is spin-j, has mass m x in the 
range 20 to 500 GeV, and is its own antiparticle. 

At present, these dark matter models are speculative. Proof will come only from direct detection of dark matter 
particles in the laboratory. Nonetheless, the models are sufficiently compelling to merit detailed examination of their 
astrophysical consequences. The remainder of this section focuses on WIMPS as they are better studied than axions. 

WIMP dark matter may be effectively collisionless today, but in the early universe WIMPs scattered rapidly with 
fermions in the relativistic plasma. Inelastic processes such as X°X° LL, where L is a lepton, maintained the 
abundance of WIMPs in chemical equilibrium until the reaction rates fell below the Hubble expansion rate (chemical 
decoupling or freezeout). This happened after the WIMPs became nonrelativistic and their abundance decreased 
exponentially through annihilation. At a temperature ~ m x /25, annihilation ceased to be effective and the abundance 
of WIMPs froze out to n(x°)/n(L) ~ 10" 10 . 

After annihilation ceased, WIMPS continued to undergo elastic scattering j°L <-* %°L with abundant leptons until 
the temperature had fallen to a few MeV (not coincidentally, slightly before neutrino decoupling at 2 MeV). During 
this period the WIMPs underwent Brownian motion as they were bombarded by the light relativistic leptons, with 
two important consequences. First, the WIMPs were maintained in kinetic (thermal) equilibrium with the relativistic 
plasma until the elastic scattering rate fell below the Hubble expansion rate (kinetic decoupling). Second, friction 
between the WIMPs and leptons damped density perturbations generated in the early universe. WIMP dark matter 
should therefore have a thermal cutoff in its perturbation spectrum at length scales corresponding roughly to the 
Hubble distance at kinetic decoupling. 

The Brownian motion of WIMPs can be studied starting from the relativistic Boltzmann equation. Expanding the 
Boltzmann collision integral in powers of the momentum transfer divided by the WIMP mass, for nonrelativistic 
WIMPs one obtains the relativistic Fokker-Planck equation 

(18) 



dx \ P »d? +r * PxPv Wi 



dx) c r dp 



(p-m x \ L )f + m x kT L j^ 



where /(x,p,f) is the distribution function in a local Lorentz frame, vl and Tl are the mean fluid velocity and 
temperature of the lepton fluid, and 7 is the collision rate coefficient determined from particle physics. 

Equation (1181 and the Einstein and fluid equations for the photon-lepton plasma have been solved numerically 
including the acoustic oscillations of the plasma before and during kinetic decoupling, the frictional damping occur- 
ring during kinetic decoupling, and the free-streaming damping occurring afterwards and throughout the radiation- 
dominated era 1 30]. For a m x = 100 GeV WIMP with bino-type interactions, kinetic decoupling occurs at a tempera- 
ture T c i = 16 MeV. Figure|2]shows the net effect of WIMP-lepton scattering on the linear transfer function of cold dark 
matter density perturbations. The damped oscillations are analogous to the acoustic peaks of the cosmic microwave 
background anisotropy and the baryon acoustic oscillations imprinted on the galaxy distribution. However, because 
WIMPs decouple at a redshift z ~ 10 10 instead of z = 1300, the length scale of these fluctuations is about 1 pc in- 
stead of 100 Mpc. This length scale corresponds to a mass 10~ 4 (r rf /10 MeV)" 3 M . It is possible, though at the 
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FIGURE 2. Cold dark matter density transfer function versus wavenumber (in units of the inverse decoupling time) after kinetic 
decoupling in the early universe; v is the factor by which primordial fluctuations grow due to gravitational instability. Here is 
the temperature at which kinetic decoupling takes place; realistic models have 0.01 < < 0.02 (larger values of this 

parameter give more damping). The upper, monotonic curve assumes, unrealistically, that the dark matter is always collisionless 
and was never coupled to the radiation. The other non-oscillating curve shows a crude model of kinetic decoupling described by a 
Gaussian cutoff. 



moment controversial, that substructures as small as this might persist in our galaxy today and be detectable through 
the products of rare dark matter annihilations 13 ill . 



BROWNIAN MOTION IN GALAXY FORMATION 

After kinetic decoupling in the early universe, dark matter dynamics is governed entirely by gravity. As shown above, 
stars undergo Brownian motion in a self-gravitating cluster. However, dark matter particles are much lighter than stars 
and their two-body interaction time is far longer than the age of the universe. Nonetheless, dark matter particles can 
scatter from quasiparticles or substructure in galaxy halos much as electrons scatter from substructure in the nucleon 
during deep inelastic scattering. In effect, galaxy halos are filled with gravitational partons. 

The concept of dark matter substructure is based on the fact that nonlinear gravitational instability accumulates mass 
into dense clumps conventionally called dark matter halos (or minihalos or microhalos for the very small ones) 1 331. In 
the hierarchical clustering paradigm for structure formation, small-scale structures form first and then aggregate into 
larger objects where they persist for some time. Gravitational N-body simulations have shown that the spherically- 
averaged radial density profiles of the evolved halos take a nearly universal form 13311 without providing an explanation 
for this result. One possibility is that scattering by substructure leads to relaxation. 

The BBGKY hierarchy provides the exact statistical description for the evolution of a classical gas. The first 
BBGKY equation is similar to the Boltzmann equation except that it is exact, not phenomenological, and it is not 
closed: 

^ +v .^ +g .|/ = _|_.F v , i-.F v =47TGp 2c (x,v,0. ( i9) 
dt dx d\ d\ dx 

Here p2 C is an integral over the phase-space two-point correlation function. The second BBGKY equation gives the 
evolution of the two-point correlation in terms of the three-point correlation, etc. For weakly correlated gases, these 
higher-order correlations may be neglected, and the two-point correlation term may be approximated by a Boltzmann 
collision integral. A different approach is needed for a gravitational plasma, where Boltzmann's Stosszahlansatz does 
not apply. 

In the early stages of gravitational clustering, p2 C can be evaluated using second-order cosmological perturbation 
theory and the first BBGKY equation reduces to a Fokker-Planck equation 1 34], with surprising results. First, the 



diffusivity eigenvalues can be negative, which appears to be a consequence of gravitational instability. Second, in 
the quasilinear regime there is no drag — y = in equation \\5\ — but there is a nonzero radial drift A = A(r)e r . 
Dynamical friction arises only in higher orders of perturbation theory (or in the fully nonlinear regime). The radial 
drift is induced by substructure. Finally, the timescale for relaxation obtained from the drift and diffusivity is the 
Hubble time, i.e. the collapse time of the initial perturbation. This initial relaxation is surprisingly fast. 

The only approximation made in deriving this Fokker-Planck equation for dark matter evolution was second-order 
perturbation theory. The starting point was the exact BBGKY equation, not the phenomenological master equation, 
and there was no assumption of Markov dynamics. Indeed, the time evolution of a given realization of the random 
process is completely smooth. Why, then, is the system described by a Fokker-Planck equation? The answer is that the 
equation of motion for dark matter particles is a modified Langevin equation: 

^ = -2H(t)v + b(t)g, (20) 
dt 

where H(t) is the Hubble expansion rate, b(t) is proportional to the growth factor of density perturbations, and g(x) is 
a Gaussian random field. For one realization of this process (i.e., one universe) this field is definite and the acceleration 
of every particle is smooth. The kinetic equation governs the evolution of the average dark matter halo, hence involves 
averaging over an ensemble of Gaussian random fields. This averaging leads to a Fokker-Planck equation. 

The derivation suggests that spatial fluctuations of the density field can cause relaxation, but is not conclusive be- 
cause the most important dynamical effects occur only in the fully nonlinear regime. Dynamical friction of substructure 
(both as "field" and "test" particles) and tidal stripping of substructure must be incorporated into the description. A 
major stumbling block is the statistical characterization of the substructure and its effects. Most investigations of sub- 
structure neglect velocity information, giving an inadequate description of phase space correlations. It remains to be 
seen whether a satisfactory kinetic theory can be devised for the fully nonlinear regime of dark matter gravitational 
clustering. 



CONCLUSIONS 

Brownian motion continues to serve as a paradigm for stochastic processes in physics and other disciplines. Einstein's 
great insight was that the dynamics of Brownian motion can be understood by applying thermodynamic equilibrium 
to the interaction between a fluid of Brownian particles and the fluid in which those particles are suspended. 

The universality of thermodynamics arises because most statistical systems near equilibrium relax at rates calculable 
from thermal equilibrium. Self-gravitating systems like globular star clusters and galaxies, with their negative specific 
heats and lack of thermodynamic equilibrium, are an important exception. While N-body simulations are the main tool 
for studying the evolution of self-gravitating systems, analytical insight can be obtained from kinetic theory, especially 
from the diffusive evolution driven by fluctuations. 

A century after Einstein's analysis of Brownian motion, the kinetic theory of self-gravitating systems remains a 
largely unsolved problem, presenting a great opportunity for a future Einstein — or a Fokker, Kramers, or Dupree. 
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